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Large N QCD - Continuum reduction 
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Numerical evidence combined with Eguchi-Kawai reduction indicate that there are no finite volumes effects 
in the large N limit of QCD as long as the linear extent of the four-torus is bigger than a critical size. This is 
referred to as continuum reduction. Since fermions in the fundamental representation are naturally quenched in 
the large N limit, as long as we only have a finite number of flavors, continuum reduction provides us with the 
exciting possibility to numerically solve large N QCD using chiral fermions and present day computers. 
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1. Introduction 

Overlap fermions enables one to compute 
quantities in the chiral limit of lattice QCD. 
But one has to work with e(75-DM,(— to)) where 
Dwi—m) is the Wilson-Dirac operator in the su- 
percritical region and e is the sign function of the 
operator. This makes the algorithms on the lat- 
tice computationally intensive and the situation is 
further complicated by the presence of very small 
eigenvalues of j5D^{—m) P]. In spite of these 
difficulties several results have been obtained in 
the quenched approximation [3j and dynamical 
simulations using overlap fermions have been ex- 
plored '5. 

Large N QCD 5^ provides an interesting alter- 
native to dynamical simulations of SU(3) QCD 
since fermions in the fundamental representation 
are naturally quenched in the — > oo limit of 
large N QCD as long as the number of flavors are 
kept finite. It is useful to use the "double line" 
notation for the gauge propagator in large N QCD 
to denote that there are two color indices associ- 
ated with the adjoint representation of the gauge 
field. In this notation, fermions in the fundamen- 
tal representation have only one line in their prop- 
agator. Only planar diagrams survive in the limit 
of iV — > oo and the vacuum diagrams arising out 
of the fermionic part are suppressed by a factor of 

and fermions are naturally quenched [S] . Large 
N QCD has all the qualitative features of iV = 3 
QCD and the spectrum of mesons has been cal- 
culated in two dimensions One also expects 



a non-trivial spectrum of mesons in the large N 
limit of four dimensional QCD but this model has 
not yet been exactly solved. 

2. Continuum reduction 

Eguchi and Kawai 7. proved that lattice QCD 
in the limit of — > oo can be reduced to QCD on 
a single site provided the global U'^{1) symmetries 
that multiplies the link matrices by U{1) phases 
are not broken. Unfortunately, this symmetry is 
broken for d > 2 8 and one cannot reproduce 
QCD in the large N limit by working on a I'' lat- 
tice. But, the proof of Eguchi and Kawai also 
holds for an L'^ lattice where L > 1. That is to 
say, the infinite space-time lattice can be reduced 
to a finite L'^ lattice in the limit of ^ cxo pro- 
vided the global U'^{1) symmetries that multiplies 
Polyakov loops by U{1) phases are not broken on 
the i'' lattice. 

The only parameter in the lattice theory is the 
't Hooft gauge coupling b = piv" This coupling 
has the dimensions of length in three dimensions 
and it is dimensionless in four dimensions. One 
can explicitly compute the location of the transi- 
tion point bc{L) such that the U'^{1) symmetries 
remain unbroken for b < bc{L). If bc{L) scales 
properly with L then one can define a critical size 
Ic in the continuum such that the [/''(I) symme- 
tries remain unbroken for I > l^. The argument of 
Eguchi and Kawai will hold for all L as long as we 
keep b < bc{L). There will be no dependence on 
the box size as long as the box size is bigger than 
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the critical size. This is referred to as continuum 
reduction 0. Physical results can be extracted 
by working on an L'^ lattice and keeping b just 
below hc{L). Computations need to be done on 
two or three different values of N to study the in- 
finite N limit. Results for two or three different 
L values can be used to study the effect of finite 
lattice spacing. 

3. The bulk transition 

One lattice artifact has to be taken into account 
when setting the parameters {L^b < bc{L)} for 
physical computation. There is a bulk transition 
on the lattice in the large N limit that is asso- 
ciated with the spectrum of the single plaquette 
(1x1 Wilson loop). There is no phase transition 
for finite N but there is a cross-over that gets 
stronger with increasing N. The eigenvalues e'^ 
of Up, the parallel transporter along the plaque- 
tte, can be anywhere in the range — tt < < it. 
But the field strength is small close to the contin- 
uum limit and is expected to be in a small range 
close to zero. In fact, there is a value of coupling 
bs such that the distribution of 9 develops a gap 
around ±7r for b > bs- This value of coupling 
does not depend on L and therefore the transi- 
tion is purely a lattice artifact. But one has to 
set b > bB to obtain the proper continuum limit. 
Only for some L > Lb does bc{L) become bigger 
than bB- For L < all one has is the bulk tran- 
sition and the [/''(I) symmetries also break at the 
same coupling. This is the reason why one cannot 
obtain continuum QCD in the large N limit us- 
ing the single site Eguchi-Kawai model and there 
is a smallest lattice size that one has to use in 
{L > LB,bB < b < bc{L)}. 

4. Numerical results 

In order to study QCD in the large N limit, one 
has to first obtain the appropriate set of parame- 
ters {L > LB,bB <b < bc{L)}. The U'^{1) sym- 
metries are not broken in d = 2 and bc{L) = oo 
for all L > 1. The bulk transition occurs at 
bB — 0.5 and Eguchi-Kawai reduction holds 
on a 1^ lattice. 

Numerical study involves the generation of 



thermalized gauge field configurations and this is 
obtained using a combination of heat-bath and 
over-relaxation updates. The heat-bath updates 
are performed on all SU(2) subgroups and the 
over-relaxation updates are performed on the full 
SU(N) group jlS . A preliminary investigation 
in three dimensions show that bs = 0.4 and 
Lb = 2 jn|. Therefore, one can obtain results 
in three dimensional large N QCD by working 
on as small a lattice as 3^. The breaking of the 
U^{1) symmetries occurs as follows. There is a 
bc{L) where one of the three C/^(l) are broken 
and the other two remain unbroken. The eigen- 
value distribution of the Polyakov lines in the un- 
broken direction are flat (random) but there is a 
peak in the distribution of the broken line. As 
one goes above bc{L) more lines gets ordered and 
there are two more transition points at any finite 
L. More investigation has to performed to study 
the behavior of the higher transitions as a func- 
tion of L, but the lowest transition corresponding 
to one ordered Polyakov line scales properly with 
L. It is quite likely that all these transitions fuse 
as one goes to the continuum limit. An analysis 
of this transition shows that 0.6 < &c(3) < 0.7, 
0.8 < &c(4) < 0.9, 1.0 < bc{5) < 1.2 and 
1.2 < bc{6) < 1.35. When combined, these re- 
sults imply 

4.2 < < 5 for d ^ 3. (1) 

Oc(L) 

A study of the twisted Eguchi-Kawai model in 
four dimensions shows that bB = 0.36 but that 
this transition is strongly first order ^21- Nu- 
merical investigation of the U'^{1) breaking tran- 
sition shows that only for L > 7 does bc{L) go 
above 0.36 ^Sl- Like in three dimensions, there 
are many transitions associated with ?7^(1) break- 
ing on the lattice. Let Oc, Ic, 2c, 3c and 4c denote 
phases with 0,1,2,3 and 4 ordered Polyakov loops 
with b above bB- Let Oh denote the phase below 
bB- One can start with a Ic phase and bring it 
down to a Oc phase and still be below b = 0.36. 
This is a metastable state. This metastable tran- 
sition from Ic Oc is shown in Figs. I1I8I The 
system is in a Ic state at & = 0.359 on a 7^ lat- 
tice with = 37 and it goes down to a Oc state 
at 6 = 0.355. We know that both b = 0.359 
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and b = 0.355 are on the weak coupling side of 
the bulk transition since the eigenvalue distribu- 
tions of the 1x1 Wilson loop shown in Fig. ^ 
and Fig. |21 exhibit strong gaps at both values 
of b. The breaking of the U'^{1) symmetry can 
be probed by looking at Polyakov loop observ- 
ables. The Polyakov loops (as well as various 
Wilson loops) were built out of "smeared" U^{x) 
matrices, rather than the original link matrices 
Ufj.{x) J3] since it suppresses ultraviolet fluctu- 
ations. Let Pfj_ be the smeared Polyakov loop 
in the /i direction and let its eigenvalues be e'^*; 
i — 1, • • • , n. The quantity |S] 

1 ^ 1 

averaged over the 3-plane perpendicular to /i will 
be close to 0.5 if the U{1) symmetry in that direc- 
tion is not broken. Breaking of C/(l) symmetry in 
the /X direction will open a gap in the spectrum 
of are p{Pfj) will drop below 0.5. A time his- 
tory of p{Pfj.) in all four directions at 6 = 0.359 
(Fig. shows that the U{1) symmetry is broken 
in one of the four directions but the time history 
at & 0.355 (Fig. Ej) shows that aU four C/(l) 
symmetries are preserved. Therefore b — 0.359 is 
in the Ic phase and b — 0.355 is in the Oc phase. 
Eigenvalue distributions of at & = 0.359 also 
exhibit the breaking of the U(l) symmetry in one 
of the directions as seen in Fig. [S] But all four 
eigenvalue distributions of at 6 = 0.355 look 
uniform as seen in Fig. |B1 Further evidence for 
b = 0.359 in the Ic phase and b = 0.355 in the 
Oc phase can be seen by looking at the eigenvalue 
distribution of the Lx L wilson loop. This distri- 
bution should not have a gap and all six orienta- 
tions of the wilson loop should look the same at 
b — 0.355 but three of the six wilson loops should 
look different from the three others a,t b — 0.359. 
This is indeed the case as demonstrated in Fig. 
and Fig. |S1 The oscillations seen in the distribu- 
tions are a finite N effect. Finally, one can see that 
the system stabilizes in the correct phase even if 
it started in the wrong phase. For example, a 
system started in the Ic phase on a 6* lattice at 
b = 0.351 ends up in the correct Oc phase after a 
few sweeps of the lattice as seen in Fig. |3 Sim- 



ilarly, a system started in the Oc phase on a 7* 
lattice at b = 0.3568 ends up in the correct Ic 
phase after a few sweeps of the lattice as seen in 
Fig.Hni 

Such a metastable state can be used to reduce L 
where one realizes large N QCD. We find that one 
can go as low as P = 5 and remain in the Oc phase. 
The Ic ^ Oc transitions are located at 6c (5) = 
0.34775 ±0.00075, bc{6) = 0.35125 ±0.00075, and 
bc{7) = 0.3564 ± 0.004. For P > 7, the transition 
from unbroken P^(l) to one broken P(l) occurs 
in the stable region and 5^(8) = 0.3605 ± 0.0015, 
6c(9) = 0.364±0.001, and 6^(10) = 0.367±0.001. 

The above values for the critical coupling in 
the range L — 5 to L — 10 can be used to study 
the scaling of the critical coupling in four dimen- 
sions 13 . The inverse function Lc{b) should scale 
according to 

if it is a truly continuum phenomenon. Since the 
above result is a two-loop result, one has to use 
the the "tadpole" jj^ replacement 

b^bi^beib) e(6) = l(rrP^,,(x)) (4) 

r / 11 \ ^ 2i^^bj 

to study scaling. Fig.llllshows that the data does 
scale properly with Pq between 0.245 and 0.275. 

5. Conclusions and future work 

The breaking of the P'^(l) symmetries for d > 2 
implies that Eguchi-Kawai reduction does not 
work on a lattice. Quenched Eguchi-Kawai 
models |H] and twisted Eguchi-Kawai models 
provide a way out of this problem. The problem 
of including fermions in the twisted version of the 
model has not been completely solved. It is pos- 
sible to extract physical quantities in the large N 
limit of QCD using the quenched version of the 
model as shown in two dimensions |16| . But it 
is necessary to sample all the N'^ minima in the 
quenched model and this could prove to be nu- 
merically difficult in four dimensions. An alterna- 
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tive to the quenched Eguchi-Kawai model is pro- 
vided by the idea of continuum reduction. Con- 
tinuum reduction imphes that there are no finite 
volume effects in the large N limit of QCD as long 
as one works with L > Lb and Bb < b < bc{L). 

The eigenvalue distribution of the single pla- 
quette exhibits a gap and the space of gauge field 
configurations naturally splits into several discon- 
nected pieces. It is natural to assume that the 
disconnected pieces correspond to different gauge 
field topologies and this can be explicitly verified 
by a numerical simulation. Since there are no fi- 
nite volume effects, one should see evidence for 
a chiral condensate if the box size is bigger than 
the critical size and the chiral condensate should 
occur due to the lowest eigenvalues scaling like 
1/N as N gets large. This also can be verified us- 
ing numerical simulations on fairly small lattices 
(L = 8-10, = 20-30). 

Since there are no finite volume effects it should 
be possible to compute meson propagators for ar- 
bitrarily small momenta as long as the box size is 
bigger than the critical size. Normally, it is not 
possible to work with momenta that are smaller 
than the inverse of the box size. This can be 
avoided in the large N limit of QCD since a con- 
stant U{1) phase multiplying all the link matrices 
in a given direction can be interpreted as mo- 
mentum carried by a quark line 8 13 il7j . This 
interpretation works as long as the U'^{1) symme- 
tries are not broken since the pure gauge sector 
is invariant under a multiplication by a constant 
phase but the fermion propagator picks up a mo- 
mentum equal to this constant phase. This is 
indeed the case if the box size is bigger than the 
critical size. A direct computation of the meson 
propagators in momentum space will lend credi- 
bility to the above argument. 
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Figure 1. Eigenvalue distribution of the 1x1 
Wilson loop at 6 = 0.359 shows a gap. 



Figure 2. Eigenvalue distribution of the 1x1 
Wilson loop at b — 0.355 shows a gap. 
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History of Polyakov loops 
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Figure 3. Time history of p{Pn) shows that the 
U{1) symmetry is broken in one of the four direc- 
tions. 
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Figure 4. Time history of p{Pfj,) shows that the 
U{1) symmetry is preserved in all four directions. 

Polyakov loop distribution 
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Figure 6. Eigenvalue distribution of at = 
0.355 shows that the U{1) symmetry is preserved 
in all four directions. 
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Figure 5. Eigenvalue distribution of at 6 = 
0.359 shows that the U{1) symmetry is broken in 
one of four directions. 



Figure 7. Eigenvalue distribution of the six dif- 
ferent orientations of the L x L wilson loop at 
b — 0.359 showing that the U{1) symmetry is 
broken in one of four directions. 
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Distribution of LXL loop 

D=4, L=7, b=0.355 
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Figure 8. Eigenvalue distribution of the six dif- 
ferent orientations of the L x L wilson loop at 
b = 0.355 indicating that the U{1) symmetry is 
preserved in all four directions. 
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Figure 10. Time history of p{Pf_i) for each direc- 
tion showing the evolution from Oc to Ic. 
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Figure 9. Time history of p(P^) for each direction Figure 11. Demonstration of scaling in fom di- 

showing the evolution from Ic to Oc. mensional large N QCD 



